m=0.8; p=3; N=2; alpha=-10;
f=@(t,Y)[-m*p/(m*(p-1)-1).*Y(1)-abs(Y(2))^((2-p)/(p-1)).*Y(2);-(N+p/(m*(p-1)-1)).*Y(2)+alpha.*Y(1).*abs(Y(1))^(1/m-1)-(1/m).*abs(Y(1))^(1/m-1).*(abs(Y(2)))^((2-p)/(p-1)).*Y(2)];
y1=linspace(-0.06,0.06,20);
y2=linspace(-0.06,0.06,20);
 [x,y]=meshgrid(y1,y2);
 u=zeros(size(x));
 v=zeros(size(x));
t=0;
for i=1:numel(x)
Yprime=f(t,[x(i);y(i)]);
u(i)=Yprime(1);
v(i)=Yprime(2);
end
quiver(x,y,u,v,'r');figure(gcf)
xlabel('y_1')
axis tight square;
hold on

[tl,yl]=ode45(f,[0,500],[-0.03;-0.06]);
plot(yl(:,1),yl(:,2),'k','LineWidth',2)

[tn,yn]=ode45(f,[0,500],[0.02;-0.06]);
plot(yn(:,1),yn(:,2),'k','LineWidth',2)
[tj,yj]=ode45(f,[0,500],[0.001;-0.06]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)

[tj,yj]=ode45(f,[0,500],[0.06;0.06]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)
[tj,yj]=ode45(f,[0,500],[0.06;0.05]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)
[tj,yj]=ode45(f,[0,500],[0.06;0.04]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)
[tj,yj]=ode45(f,[0,500],[0.06;0.03]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)
[tj,yj]=ode45(f,[0,500],[0.048;0.06]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)
[tj,yj]=ode45(f,[0,500],[0.049;0.06]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)
[tj,yj]=ode45(f,[0,500],[-0.049;-0.06]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)
[tj,yj]=ode45(f,[0,500],[0.035;0.06]);
plot(yj(:,1),yj(:,2),'k','LineWidth',2)
